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Abstract 

On the basis of Langevin equation the optimal SUSY field scheme is formulated to discribe 
a non-equilibrium thermodynamic system with quenched disorder and non-ergodicity effects. 
Thermodynamic and isothermal susceptibilities, memory parameter and irreversible response 
are determined at different temperatures and quenched disorder intensities. 
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1 Introduction 

Recently, the microscopic theory of non-equilibrium thermodynamic systems with bro- 
ken ergodicity and exhibiting the memory effect, has been the subject of major interest. 
Spin glasses and random heteropolymers @, that have received much consideration, 
are well-known examples of such systems. Despite the bulk of theoretical studies had 
been employed the replica method to approach the problem analytically, there is increas- 
ing interest in alternative methods that go beyond the replica trick. The supersymmetry 
(SUSY) approach, evolved within stochastic dynamics theory governed by Langevin equa- 
tion, gives a good example of method of such kind. 

According to this method, generating functional of Langevin dynamical system is rep- 
resented as a functional integral over the superfields with Euclidean action by means of 
introducing Grassmann anticommutating variables. These variables and their products 
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serve as a basis for superfields with components that involve Grassmann fields along with 
real(complex)-valued ones. 

As it was shown in ||, static replica treatment of spin systems bears striking similarity 
with the dynamics expressed in terms of superspace within the framework of the SUSY 
method. The latter is based on using of nilpotent variables. Two-point correlator of SUSY 
field can be written in the form of expansion with coefficients that give correlators of 
observables such as structure factor S, and retarded and advanced Green functions G T . 
The memory and nonergodicity effects are allowed for by incorporating the additional 
terms q, A into the correlators. The resulting self-consistent SUSY scheme gives a set of 
equations for memory q and nonergodicity A parameters to be determined as functions 
of temperature T and quenched disorder h. 

For SUSY scheme formulation SUSY as a gauge field needs to be reduced to irreducible 
components. By analogy with electromagnetic field, that can be splitted into vector and 
scalar fields, 4-component SUSY field can be divided into chiral components that con- 
sist of regular and Grassmann constituents |3J]. In Sections 2 and 3 SUSY field will be 
reduced to 2-component nilpotent field. The latter has an advantage over conventional 
SUSY representation because its components have an explicit physical meaning of order 
parameter and conjugate field (or amplitude of its fluctuation). This rises the question 
as to optimal choice of the basis for making expansion of SUSY correlators. Currently, 
two types of such basis are known [§,0. The first one contains 3 components: advanced 
and retarded Green functions G± and structure factor S. The second basis corresponds 
to the proper 4-component SUSY field and contains 5 components which, in addition 
to the above mentioned ones, include a couple of mutually conjugated correlators of the 
Grassmann fields. In Sect. 4 it will be shown that the second basis can be reduced to the 
first one. 

The work is organized as follows. In Sect. 2, the simplest field scheme is formulated in 
terms of the 2-component nilpotent fields with second component is taken to be either an 
amplitude of fluctuation or a conjugate force (see subsections 2.1 and 2.2). In Sect. 3 the 
above-discussed method for reduction of the 4-component proper SUSY field to different 
2-component forms is presented. In Sect. 4 we show that the reduction results in decrease of 
the number of components of the SUSY correlator basis, due to the fact that the conjugate 
correlators of the Grassmann fields are equal to the retarded Green function in accordance 
with Ward identities. The SUSY perturbation theory, stated for both cubic and quadratic 
anharmonicities in Sect. 5, makes expressions for the SUSY self-energy function simple 
to calculate. This function enters the SUSY Dyson equation derived in Sect. 6 on the 
basis of effective Lagrangians for both thermodynamic systems with quenched disorder 
and random heteropolymers. Nonergodicity and memory effects are investigated in Sect. 7. 
The corresponding self-consistent equations are obtained. Behavior of a non-equilibrium 
thermodynamic system for various values of temperature and quenched disorder intensity 
is analyzed in Sect. 8. Appendices A, B, C provide details concerning the SUSY formalism 
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under consideration. 



2 Two— component SUSY representation 



Let us start with the simplest stochastic Langevin equation governing the spatiotem- 
poral evolution of order parameter r)(r,t): 

r)(r,t) - DV 2 V = -<y(dV/drj) +CM), (1) 

where the dot stands for derivative with respect to time, V = d/dr, D is the diffusion- 
like coefficient, 7 is the kinetic coefficient, V{rj) is the synergetic potential (Landau free 
energy), C( r ;^0 is a Gaussian stochastic function subjected to the white noise conditions 

(C(r, t)) = 0, (C(r, t)C(0, 0)) = ^T5(v)5(t), (2) 

where the angular brackets with subscript denote averaging over the Gaussian proba- 
bility distribution of (, T is the intensity of the noise (the temperature of thermostat). 

Further, it is convenient to measure time t, coordinate r, synergetic potential V, and 
stochastic variable C, in units t s = ( 7 T) 2 /D 3 , r s = jT/D, V s = D 3 /7 3 T 2 , ( s = D 3 /(-fT) 2 
respectively. The equation of motion (]I|) then reads 

fi{r,t) = -6V/6ri + C(T,t), (3) 

where short notation is used for the variational derivative 



5V/5r] = 5V{r]}/57] = dV{rj)/dri - V 2 r?, V{r]} = J 



V(V) + \ (Vr^) 2 



dr, (4) 



the coefficient 7T in Eq.(|[) becomes unity and the distribution of variable ( takes the 
Gaussian form 

P {C}ocexp(-i|c 2 (r,t)drdt). (5) 
The basis for construction of the field scheme is the generating functional || 

Z{u(r,t)} = Jz{r]}exp (^Jurjdrdt^Dr], (6) 



^(r,t)} = (|n 5 {^ + ^-c}det 



6r) 



(7) 



so that its variational derivatives with respect to an auxiliary field u(r,t) give correlators 
of observables (see Eq. (|T2])). Obviously, Z{u} represents the functional Laplace transfor- 
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mation of the dependence Z{r]}, 5-function reflects the condition (|3|), the determinant is 
Jacobian of the integration variable change from ( to rj. 



2. 1 Fluctuation amplitude as a component of nilpotent field 



Further development of the field scheme proceed depending on the type of connection 
between stochastic variable ( and order parameter 7]. For thermodynamic system, where 
the thermostat state does not depend on rj, the determinant in Eq.(|7]) assumes constant 
value that can be chosen as unity. Then, by using integral representation for 5-function 



lOO 

§{x(r, t)} = J exp ^— J Lpxdrdt^J Dip 



with the ghost field (p (r, t) and averaging over distribution fl5|), we have the functional (|7|) 
in the standard form 

Z{ V (r,t)} = J exp[S{7i(T,t)Mr,t)}]D<p, (9) 

where the action S = J £drdt is measured in units S s = 7 2 (T/D) 3 with the Lagrangian 
given by 

L( v , tp) = (iprj - V 2 /2) + <p{5V/6r,). (10) 

In order to obtain a canonical form of the Lagrangian ( |H1| ) let us introduce the nilpotent 
field 

(f> v = -n + ti(p (11) 
with Bose components rj, tp, and nilpotent coordinate $ obeys the relations 

■Q 2 = 0, Jd$ = } J tfdtf = 1. (12) 

As is shown in Appendix A, the first bracketed expression in Lagrangian ( ID ) takes the 
form of kinetic energy in the Dirac field scheme : 

K = - J <j)D(j)d§. (13a) 
Hereafter indexes are suppressed. The Hermite operator D is defined by equality 

D ' = -^ + ( 1 - a »Si)i (14) 
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and enjoys the property ( |A.6| ). On the other hand, the algebraic properties (0) of coordi- 



nate § allow to rewrite the last term in Eq.QlOl) in the standard form of potential energy 
(see Appendix A) 

7T = j V{(j))d!&. (13b) 
The resulting Lagrangian ( |T0D of the Euclidean field theory is 

L = k + 7t = J X di?, \(<f>) = ^4>D(f) + V(4>). (15) 



According to Appendix A, the expressions (P~D[) , (|T3p become invariant with respect to 
transformation e eD given by operator ([14]) if only a parameter e — > is pure imaginary 
and the fields 77(1*, t), <p(r, t) are complex-valued. Then, operator D is the generator of the 
nilpotent group. 

After equating the first variation of the functional 

S {0(C)} = [ A(0(C))dC, ( = {r,t,$} (16) 



to zero, we obtain the Euler-Lagrange equation 
5X 5X 

D m + rr a ' (17) 

Substituting the expression flT5| ) in Eq.(|T^) yields the equation of motion 

D(f) + 5V/5(j) = 0. (18) 

Projection along axes of usual and nilpotent variables gives the system of the equations 
f] = -SV/Sr] + (p, (19) 
6 2 V 

= ¥>, (20) 
o z r] 

that determines kinetics of the phase transition. Being obtained from the extremum condi- 
tion for Lagrangian (|i~0|) these equations determine the maximum value of the probability 
distribution 

P{r / (r,t),^(r,t)} = Z- 1 exp(- J L( V , ^>)drdt) , (21) 

that specifies the partition function Z = Z{u = 0} in Eq.(^|). Comparison of expression 
(|j~9l ) with Langevin equation (^|) leads to the conclusion that the quantity if determines 
the most probable value of fluctuation of the field conjugated to the order parameter. On 
the other hand, it means that the initial one-modal distribution (|5]) transforms into the 
final two-modal form (PT|). 
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2.2 Conjugate field as a component of nilpotent field 



It is well to bear in mind that there is another representation of two-component nilpotent 
field. Let us introduce field f(r, t) defined by the relation 

fi = f + V>. (22) 

Then the Lagrangian fllO|) takes the form 

Since the last term of Eq.(^) is the total derivative of V with respect to time, its contri- 
bution to the partition function gives a factor that is integral over initial and final fields 
77i(r) = rj(r, ti), rjf(r) = r)(r,tf) (here we return to dimensional magnitude of the potential 
V). 

z = j ew (_vM^XMj DviDvf (24) 

The remaining part of Lagrangian fl23|) yields the Euler equations 
5 2 V 

1 = -J^f, <*) 
/ = -8V/5r}. (26) 
Differentiating Eq.(|19|) with respect to time and taking into account Eqs.(pOl), fl22|), it is 



not difficult to derive Eq.(|25|). As for Eq. (|26l ), it defines f(r,t) as the field conjugated 
to the order parameter r](r,t). Note that Eq.(p6"D implies the force / explicitly does not 
depend on the time t. 



By analogy with the definition (|TTJ) let us introduce now another nilpotent field || 

<P f =V-$f, (27) 

where Bose components are the order parameter rj and the force / with opposite sign. As 
it is shown in Appendix A, expression for Lagrangian in terms of 0/ has the same form 
as in Eq . (|15D with the generator of the nilpotent group given by 

Note that Df obeys the algebraic relation (|A.6|). 
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2.3 Connection between two -component nilpotent representations 



In this subsection we discuss the relation between the two above two-component nilpotent 
fields (|TT| ) and (27) that makes using of the fields algebraically equivalent. 

Let us introduce the operators t± = e ±,&dt , d t = d/dt that induce the following transfor- 
mations of the fields 4> v and 0/ 

r±(p TV (t) = T/ (t), r ± <f> ±f {t) = <f> ±<p (t). (29) 



Eq.(P9"D shows that operators t± transform the field to its counterpart. So we have the 
mappings relating the representations. 



By making expansion in power series over with help of Eqs.([12D, (|22[ ) one obtaines 
<P TV (t ± 0) = <P Tf (t), ±/ (t ± 0) = <f) ±v (t), (30) 

that shows that operators t± shift the physical time t by the nilpotent values ±$: 

T±(f>^(t) = <f>^(t ± 0), r±0 ±/ (t) = <P± f {t ± 0). (31) 

The same results can be obtained by using matrix representation defined by Eqs. ( |A.7|) , 
(PO| ) and ( PTTOQ . 

On the other hand, the above mappings r + (pf = <p v , T-0 V = 0/ induce corresponding 
transformations of the generators (|HD , ( PH|) 

D, = r_^r + , D v = r+D/T- (32) 



Note that the action with Lagrangian (|Tq) is covariant with respect to transformations 
fl29|), provided / = (potential 1/ does not depend on time explicitly). 



3 Reduction of proper SUSY fields to the two— component forms 



The considerations given in previous section rest on the assumption that the Jacobian 
of variable change from ( to the order parameter rj is constant. However, in general case 
determinant of an arbitrary matrix \A\ can be expressed as an integral over Grassmann 
conjugate fields ip(r,t), ip(r,t), that meet conditions type of Eqs. fll"2]) 

det \A\ = J exp ($Aip) d 2 ip, d 2 ^ = dipdrp. (33) 

Physically, the appearance of new degrees of freedom ip, tp means that the state of thermo- 
stat turns out to be dependent on the order parameter — as it is inherent in self-organized 
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system []10[ . As a result, the Lagrangian (|T0|) supplemented with the Grassmann fields ■?/>, 



■0 takes the form 

/ 5V \ -(d 5 2 V\ 

C( V , ^).^- T + _ v j-^- + _j^. (34) 

Introducing the four-component SUSY field 

® v = V + W + '$0 + 09(p, (35) 
by analogy with previous section the SUSY Lagrangian is 

£ = J Ad 2 9, A(%) = ^{V v %) (V^ v ) + V(%), d 2 9 = d9d9, (36) 

where 9, 9 are Grassmann conjugate coordinates that replace the nilpotent one As 
compared with Eq . (fl5|) , where the kernel A is linear in the generator flT4j ), a couple of the 
Grassmann non-conjugated operators 

D *=|r 2 4' p *=ss (37) 

enters the expression for SUSY Lagrangian. The Euler equation for SUSY action reads 
-i|^,P]<& + ^ = 0, (38) 



where the square brackets denote the commutator. Projection of Eq. (|38|) along the SUSY 
axes 1, 9, 9, 99 gives the equations of motion 

fj - V 2 i] = -dV/drj + <p, (39a) 

+ vV = (d 2 V/dT] 2 )<p - (d 3 V/dr] 3 )^ij, (39b) 

^ - vV = -(d 2 V/d7] 2 )ij } (39c) 

- V 2 ^ = -(<9V/«9r/ 2 )^, (39d) 

that give Eqs.(^), (|20|) at ^ = ^ = 0. It can be readily shown that this system can be 
obtained from the Lagrangian fl3"4|). From Eqs. ( |39cj ) and ( |39d[ ) we obtain the conservation 
law S 1 + Vj = for the quantities 

S = ^, j = (V^-^(W>). (40) 

For inhomogeneous thermodynamic systems S is a density of sharp boundaries, j is a 
corresponding current ||. In particular, the approach of the four-component SUSY field 
complies with the strong segregation limit requirement of copolymer theory []TT | . For self- 



organized system the magnitude S gives the entropy, j is the probability current JTDI . 
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So, for thermodynamic system, where the entropy is conserved, we could disregard the 
Grassmann fields ip(r,t), ip(r,t) = const. As a result, the four-component SUSY field 
( |35| ) is reduced to the two-component form (|TT|). 



In order to justify this let us write the kinetic term of the Lagrangian ( |36|) in the form 
"(1/4)** " 



V V 



§ v where 



V V 



2 



8989 



+ 1-20 



d_\ d_ 
89! di' 



(41) 



The expression ( |41| ) restricted to two-component form with d = 99 yields the generator 
([14] ) as is needed. It is of interest to note that variable d satisfies (|T2|) . In addition, since 
the self-conjugated value $ = $ is commutating quantity, it is nilpotent rather than 
Grassmannian. 



As in the case of two-component nilpotent fields in Section II, one can go over from the 



fluctuation amplitude <p to the conjugate force / by using Eq. (|22|) . Then, the first bracket 
in Lagrangian fl3~4"D takes the form (p3|) and instead of the system (39) one obtains the 
equation (cf. Eq.fl2"5|)) 

f, = -(S 2 V/5n 2 )f - {5 3 V/5n 3 )^ij (42) 

supplemented with the definition of force (|26|) and the equations ( [39c| ,d) for the Grassmann 
fields tp(r, t), ip(r, t). As above, the equation of motion (^) can be derive by differentiating 
Eq . ( |39a| ) with respect to time and taking into account Eqs.(f22"D, ( |39b|) . The corresponding 
Lagrangian C(r], f, ip, ip) takes the SUSY form (cf. Eqs.(|3II)) 



C = J Ad 2 fl, 
with the SUSY field (cf. Eq. 



(43) 



$ / = n + 9i\) + ip9 - 99 f = $ v - 66$ f = T_$ v , 

T± = e ±m \ d t = d/dt (44) 



and the Grassmann conjugated operators (cf. Eqs.(|3"7|)) 

9 8 _ d 



V f=xE-^ V f = --9-. (45) 



By analogy with Eqs. (|29D -(|3T]) with operators r± = e ±,ddt replaced by T± = e ±e99t , where 
d t = d/dt, the SUSY fields (|3~5"D, ([44]) can be transformed into each other and couples of 
the SUSY operators fl3~7|), ([45"D are related by means of transformations (cf. Eqs. (|32|) ) : 

V f = T„V V T + , V f = T_V^T + . (46) 
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According to Eq. fl45|) , kernel of kinetic part of the SUSY Lagrangian (fl"3| ) is (cf. Eq.([41|)) 



-^-(ll+^W4-fe (47 » 



Note that the operator (28) can be obtained from Eq.([|7]) by taking into account the 
condition of the Fermion number conservation 9(d/d9) = 9(d/d9) and by setting 99 equal 
to •&. 



So, both four-component Grassmann fields (p5[) , ([44]) with SUSY generators given by 
Eqs.(|3T|), (|45| ) can be reduced to the corresponding two-component fields, Eqs.flTTD, (p7|), 
with operators (|14|) , (p8|), respectively 

It is worthwhile to mention that such reduction can be obtained according to the SUSY 
gauge conditions 

£>$ = 0; P$ = 0. (48) 



Indeed, according to definitions (|35|), (137|), fl44]), (f45|) the equalities ( |48D give the relation 



#V + # - 200/ = 0, (49) 



that reduces the SUSY field ([44]) to the form (p7D with opposite sign before /, provided 
d = 99. 

Despite of the same number of components, one has to have in mind that the reduced 
SUSY field from Eq.(p7[) and couple of Grassmann conjugate chiral SUSY fields ( |B.9| ), 
which appearance is a consequence of SUSY gauge invariance also (see Appendix B), have 
different physical meaning. The main distinction is that the first field consists of two Bose 
components rj, f, whereas the chiral SUSY fields <f>+, 0_ are the combinations of Bose i] 
and Fermi ip, ip components. Formally, this is due to the fact that for separation of the 
chiral SUSY fields the conditions ( |B.7| ) of the SUSY gauge invariance are fulfilled not 
for the initial SUSY field $, which satisfies to conditions (f4"8[), but for components $±, 
resulting from $ under the action of operators T± = exp (±99d t ) (see Eq.(pp). 



According to the above considerations, the transformation operators T±, that shift physi- 
cal time t by Grassmann values ±99, relate the SUSY fields (|35|) , ([44]) and corresponding 
generators (|37|) , (|45|) . It should be emphasized that only the latter form a pair of Grass- 
mann conjugated operators. The physical reason of this symmetry is that the correspond- 
ing equation of ([42]) is invariant with respect to the time inversion, whereas the equations 



( |39a| ), ( |39b| ) for components of the SUSY field ([35]) are not. However, in addition to the 
field = $ + obtained from the initial field $/ under the action of operator T + , another 
SUSY field $_ emerge under the action of operator T_ that shifts the time t in opposite 
direction. From Eqs.( p.5| ), (|22| ) it can be seen that the fields <1>-|- = ^^(it) correspond to 



opposite directions of time. 
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However, equations (|39c| ), (|39d| ) for the Grassmann components ip(r,t), ip(r,t) are invari- 
ant under the action of T± To break the invariance let us introduce additional operators 
of transformation 

f± = exp [e (5±9ip + <5 T #)] (50) 

where source parameter e — > 0; 8 + — 1, <L. = for the positive time direction and 5 + = 0, 
5- = 1 otherwise. The Euler SUSY equation (|38|) for transformed superfield <3>± = T±^ v 
is reduced to the components 

■fj — V 2 rj = —dV/dr] + (p — etptp, (51a) 

(p + vV = (d 2 V/dr} 2 )ip - {^V/drf)^ + e^ip, (51b) 

^ - V 2 V> = -(9V/^ 2 )V - £ {8-{$/il))r) + 8+ [(fi -(p) + (d 2 V/dr] 2 )r]\ } V, (51c) 



t/j + V> = (9V/9r/ 2 )^ - e {5+(ip/ip)ri - 5. [(fj - <p) + (d 2 V/dr] 2 )r]\ | ip, (51d) 

where the terms of first order e are kept. These equations give Eqs.(39) at e —>■ 0, but 
combination of Eqs. (|51c|) , (|51d| ) at e ^ leads to the following equation for the quantities 

5 + Vj = ±eFS, F = dV/dr] - 2(d 2 V/d7] 2 )r] (52) 

instead of the law of entropy conservation. Since entropy S of a closed system (Vj = 0) 
increases in time, provided F > 0, in Eq . (|52"D one has to choose the upper sign correspond- 
ing to the positive time direction. So, the operator ( |5T^ ) breaks symmetry with respect to 
the time reversibility. The above-mentioned condition of positiveness for effective force 
F = dV/drj — 2(d 2 V/dr] 2 )r] means that the effective potential V is an increasing convex 
function of the rj that is inherent in an unstable system. It is of interest to note that near 
the equilibrium state, where dV/drj — 0, rj <C 1, the force F ~ —(d 2 V/dr] 2 )r] is always 
positive for unstable system. 



Finally, in order to visualize the difference between two-component nilpotent fields (|lT|), 
P7| ) and chiral fields ( |B.9j ) let us represent the SUSY field (H) as a vector in four- 
dimensional space with axes 6° = F = 1,6,6,96 = Then conditions © of the SUSY 
gauge invariance mean that field (|44|) is reduced to the vector ( p7|) belonging to a plane 
formed by axes 1, Accordingly, the conditions (|B.7|) of the chiral gauge invariance split 
total SUSY space into a couple of orthogonal subspaces, the first of which has the axes 1, 
6 and contains the vector 0_, and second — axes 1, 6 and vector <f> + . Since these subspaces 
are Grassmann conjugated, <ft_ = 0+, it is enough to use one of them, considering either 
vector </>_, or (f> + (see Appendix B). Such program was realized in Ref.|12j, whereas the 
above used nilpotent field (^7|) is derived by projecting chiral vectors <f>± to a plane formed 
by axes 1, It follows that our approach stated on the using nilpotent fields (fTCD, (J27T) 
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and the theory |Ll| are equivalent. The SUSY method presented in the book |13] is also 



based on usage of the chiral fields 0_ = (p — iip6, <p + = r] + 6ip (cf. with ( B.9 )) that contain 



the fluctuation ip as a Bose component of the field 0_ and the order parameter r\ in field 



4 SUSY correlation techniques 



In this section correlators of the proper SUSY fields (|35|) , fl44|) will be studied. It will be 
shown how the relevant correlation techniques can be reduced to the simplest scheme by 
making use of the two-component field (|Tl|). 

To begin with let us introduce the SUSY correlator 

C(z, z') = ($(z)$(z')>, z = {r, t, 6, 6}. (53) 

From the equation of motion (|38|) we have the equation for Fourier transform of the bare 
SUSY correlator C^°\z,z') with the potential Vq = (l/2)<3> 2 in the following form 

L^(6)C { ^(6,6') = 5(6,6'), L = 1 — (1/2)[V,V], (54) 
where 5(6, 6') is the Grassmann 5-function 

s(e,e') = {e-e){6-6'), (55) 

u is the frequency and k is the wave vector. The solution of Eq.([54]) reads 



_ y iHU2)W,v})me') 

where the indexes u, k are suppressed for brevity. From the definitions (|3"7|), ( |5^| ) and 
equality ^D,Vf = -Auo 2 (see Eqs. flPlD ), the bare SUSY correlator for SUSY field (||) 
can be written in the explicit form 

C (o) (g> n _ 1 + (i - M(g - g> - (i + j^Kg - gy (57) 

In the case of the SUSY field (|44D , by using transformation (]46[) the above result is found 



to be modified by adding the term iuj(66 — 6 6') to the numerator of Eq. (|57|) . 

It is convenient to introduce the following components as a basis for expansion of SUSY 
correlators 
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T(e, e') = 1, B (e, 9') = ee, b 1 (9, e') = e'e', (58) 

F (6, 9') = -e'e, Fi(6, 9') = -99'. 
Let us define the operator product 

X(9, 6') = J Y(0, 9")Z(9", 6')d 2 6" (59) 

for superspace functions Y, Z. Eg . (|59|) immediately provide the multiplication rules for 
the basis operators ( |58"D summarized in Table I: 



Table I 



l\r 


T 


B 


Bi 


F 


Fi 


T 





T 











B 





B 











Bi 


T 





Bi 








F 











F 





Fi 














Fi 



The operators T, B 0) i, F .i then form the closed basis, so that expansions for correlators 
are (see Eqs.flCTD, (|6lOD ) 

C v = ST + G + (B + F ) + G_ (B x + F x ) , (60) 
C f = ST + m+B + m_B! + G+F + G_Fi 

where in accordance with Ward identity (|C.6|) corresponding to the first generator (|C.5|) 
term proportional to 999 9' is dropped. Inserting SUSY fields ([35]), (0) into Eq.(|53D 
provides the coefficients of expansions (|60| ) (cf. Eqs. (|C.8| ), ( |C.11|) ): 



S = (|r/| 2 ); m + = (r]*)f ext , m_ = (r])f* xt , f ext = -/; (61) 

So, quantity S is the autocorrelator of order parameter r] and magnitudes m T meet the 
condition m* + = m_ and determine the averaged order parameter (77) corresponding to 
external force / ext = — /. The retarded and advanced Green functions G T give the response 
of order parameter 77 to fluctuation amplitude ip and vice versa (moreover, functions G± 
determine correlation of the Grassmann fields ip, ip). As it is known 0], the Fourier 
transforms G T (u) of retarded and advanced Green functions are analytical in upper and 
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lower half-planes of complex frequency u with cut along real axis u/. There is the jump 
G-(u') — G + {uj') = 4i ImG_(c</), so that the relations QC.9D , ( |C.12| ) assume the usual 
form of the fluctuation-dissipation theorem: 



G ± (u) = m±(u) =F iuS(u), S(u r ) = (2/ J) ImG_(w') 
where the frequency uj' is real. The expression for bare correlator (|57|) gives: 



(62) 



S (0) 



(0) 

m± 



(1 + u 



2\-l 



gL 0) 



(63) 



Integrate the last equation of (|62|) and taking into account the spectral representation 
~ du/lmC(u') 



C(uj) 



we arrive at useful relation 

S{t = Q)=G±(u = 0) = X , 
where the last identity is the definition of susceptibility x- 



(64) 



(65) 



The expansions (|60|) make it possible to handle the SUSY correlator (|53D as a vector of 
space constructed as the direct product of the SUSY fields (pJH) or (H3). The representation 



(|3~5p is of special convenience because it allows using of reduced basis 

A = Bo + Fo, B = Bi + Fi. (66) 
Along with T, they form more compact basis and obey the following multiplication rules: 

Table II 



l\r 


T 


A 


B 


T 





T 





A 





A 





B 


T 





B 



The expansion Eqs.(^Dj) then takes the form 



C 



ST + G + A + G B. 



(67) 



So, using Ward identities allows to get rid of autocorrelators of the Grassmann fields 
ip, ip (see relations (|5"ID). As a result, there are three basic correlators: the advanced 
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and retarded Green functions G± and structure factor S. They yield the most compact 
expansion (p7|) for arbitrary SUSY correlator of fields fl35|) . It is ready to show that 
expansion of the same form can be obtained on the basis of the two-component field QTTj) 
representation. Indeed, in this case by comparison between equations of motion (pl|) and 



58|) the commutator — (1/2) [P, V] in expression (RB) should be replaced by generator (fT3) 



and nilpotent 5-function should be 5($ — $')='& + So the resulting bare correlator is 
C'°»K^)^ 1 + (1 - Mtf+ < 1 + i ^' (68) 

1 + U) 1 

instead of Eq.(^). It is easily to see that using the definitions (cf. Eqs.(|58|)) 

r(0,0') = l, A{&,-&')=#, B(&, •&') = ■&' (69) 

gives the relevant expansion ( |S7|) . As a result, in what follows we can use two-component 
field flTTl). 



In particular, for inverse of the SUSY correlator (|67|) we have 

CT 1 = -G^SGZ^T + G^A + G^B. (70) 



It is worthwhile to note that according to definitions (pB|), (|5^D the basis operators A = 
B = provided 6 = 9', so that C(9, 9) = C(&, ■&) = S and 

J C{z, z)dz = J S{r, t; r, t)drdtd 2 9 = 0, (71) 
j C(C, C)dC = / S(r, t; r, t)drdtd$ = 

where z, ( are sets of variables {r,t,9,9}, {r,t, respectively. In the diagrammatic 
representation identities fl7ID imply the absence of bubble graphs contribution. The latter 
considerably reduces the number of graphs contributing to expansion of the perturbation 
theory (see below). 



5 SUSY Perturbation Theory 



Let us begin with the formula 

t 2 Z{u(()} 



c(C,C) 



u=0 



HOHC) _ 

where generating functional (see Eqs.©, @) 
Z{u} = ( exp ( / (pud( 



(72) 



(73) 
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has the form of average over distribution (cf. Eq. fl2T|) ) 

P{0} = Z' 1 exp (S{(f)}) , S{<j>} = J A(0)dC, (74) 

with the Lagrangian A defined by Eq. (|15D . In the zero-order approximation the action is 
quadratic 

So = \j (pLtpdC, L = l + D, (75) 

where generator D is given by Eq.([L4]). Corresponding distribution takes the SUSY Gaus- 
sian form (cf. Eq.(|D) 



1/2 



W} = (^^) ex P WdC} • (76) 

So for the bare supercorrelator we have the expression 

C(°\(,(') = L- 1 6((,C), 6(4,tf)=4 + tf, (77) 

that leads to Eq.flU) with -(1/2)[P,Z>] replaced by D if Eqs.flFD, (|TJ) are taken into 
account. The linear operator L = (C^) -1 in accordance with Eq.(|7(]) takes the form: 

L = LT + L + A + L_B; (78) 
L = —1, L± = 1 ± iuo. 



To proceed, one need to separate out anharmonic part of exponent in distribution 

fl73| ) as a perturbation and to make expansion in power series over S\. Insertion of this 
series into Eq.(|72"D gives 



C(C, CO = E MO 0W» n #0>„ > ( 79 ) 

n=0 U - 

where subscript "0" means averaging over the bare distribution (|75|). Further one has to 
make factorization by making use of the Wick theorem. Then within the n-th order of 
perturbation theory the expression ( |79| ) takes the form 

c^(CC') = JI C^(cm {n) (CiX2)C^(C2,C')dCi dC 2 , (so) 

where EW(Ci, C2) is the SUSY self-energy function of n-th order that should be calculated. 
The result essentially depends on the form of Vx(<f>) that describes self-action effects. In 
what follows we will analyse two widely used models. 
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5.1 4 -model 

Let the self-action potential be defined by the quartic dependence 

^i(C) = ^0 4 (C), C = {r,t,tf} (81) 

with the anharmonicity constant A > 0. Then terms of the first and second orders of series 
© are 

C {1) (CC) = ~j (0(C)(0(Ci))V(C')) o dCi, (82a) 

C {2) {CX') = j { (~) I / (0(C)(0(Ci)) 4 (0(C 2 ))V(C')) o dCi dC 2 . (82b) 

Now one has to count the number of possible pairings when using the Wick theorem. In 
Eq.( |S2a| ) the total number of pairings is 12, and the formula (|82a|) reads 

cP-HCO = -\j c^(c,Ci)c i0 \Ci,Ci)c^(Ci,C)dCi = o, (83) 

where Eqs.(|53|), ( |7l~D are taken into account. In Eq. ( |82b| ) the total number of pairings 
equals 192 and the Wick theorem gives 

C (2) (C,C) = y // C<°)(C,Ci) (cW(Ci,C 2 )) 3 C(°)(C 2 ,C / )dCi dC 2 . (84) 

Then, in accord with Eq . (|80D the SUSY self-energy function in the second order of per- 
turbation theory reads 

^CC'H^^CC')) 3 - (85) 

o 

Here in terms of usual diagram ideology bare correlator is replaced by exact one. 

In the diagrammatic representation terms (82a,b) correspond to the following graphs: 



According to the rule (|7l|) the former does not contribute to correlator, whereas the latter 
does 
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By analogy with the SUSY correlator fl67|) it is convenient to expand the SUSY self- 
energy: 



£T + E+A + £_B. 



(86) 



To determine the coefficients E±, £ it should be taken into account that the multiplication 
rules in Eq.(^) differ from ones given by Table II. The reason is that Eq.(^) contains 
"element-to-element" products of nilpotent quantities instead of the above operator 
product. Hence one has to use the multiplication rules given by the Table III: 

Table III 



l\r 












A(tf,tf') 


£(#,#') 




A(tf,tf') 










B (-&,-&') 









As a result, the coefficients of expansion (|86| ) take the form: 
E(f) = (A 2 /6)5 3 (t), 
E±(f) = (A 2 /2)S 2 (t)G ± (t). 

In the frequency representation that will be needed below we have 
A 2 f du)i du) 2 



E+ M = — 



(3 y (2tt) 2 
A 2 /" dcui dco>2 



S(cj — Ui — lu 2 )S(ui)S(lu 2 ), 
G±(uj -uj x - u} 2 )S(ui)S(u} 2 ). 



(87a) 
(87b) 



(88a) 
(88b) 



The obvious inconvenience of this expressions is the presence of convolutions. To get rid 
of them let us use the fluctuation-dissipation theorem 



E(t = 0) = S ± (cj = 0) 
in the form of Eq.(^). Then from Eqs. ( |S7a 
S±(w = 0) = (A 2 /6) X 3 . 



one obtains: 



J9) 



(90a) 
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5.2 Cubic anharmonicity 



Apart from the 4 -model studied above, there is a number of physical systems type of 
copolymers || where cubic anharmonicity 

^(0 = ^(0, C = {r,t,$} (91) 

has a dominant role ( fi is the anharmonicity parameter). By analogy with Eqs.(82) it can 
be shown that the first non-vanishing contribution to the SUSY correlator (jSOD is 

C {2 \C,C) = ^ ("1)7/ <^(O(0(Ci)) 3 (0(C 2 )) 3 0(C')) o dCi dC 2 - (92) 

To facilitate the factorization of these products let us depict possible graphs of the second 
order in cubic anharmonicity \x: 



The first of these graphs contains the bubble and does not contribute to the correlator. 
The contribution of the second graph is 



2 



c<°)(C,Ci) (c (0) (Ci,C2)) 2 c(°)(C 2 ,C / )dCi d( 2 . (93) 

As a result, the SUSY self-energy function reads 



E(C,C0 = y(C(C,0) a . ( 94 ) 

where the bare SUSY correlators are replaced by exact ones. By using the multiplication 
rules from Table III the coefficients of the expansion (|8SD are derived 



£(£) = Gu72)S 2 (t), (95a) 
£±(t) = (i 2 S{t)G±{t). (95b) 



These expressions, combined with Eqs.(87), determine the SUSY self-energy function 
completely. By analogy with Eq.(|90a|) we have the relation 



£ ± (cj = 0) = £(t = 0) = (fi 2 /2) X 2 (90b) 
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Finally, the resulting expressions for coefficients of expansion ( [36|) with both cubic and 
quartic anharmonicities included are 

£W = ^ 2 + yS(t))s 2 (t), (96a) 

S±(t) = (V + ^S(t)j S(t)G ± (t), (96b) 

S ± (^ = 0) = V/i 2 + yx)x 2 (96c) 



6 Self— consistent approach 

6.1 Effective SUSY Lagrangian 

Let us start with the total SUSY action taken in the site representation: 

S = S + S 1 + S int ; (97) 

S o = J 0zM)[l+£(tf)]0zM)dtdtf, (97a) 

St = J2f V 1 (<f>i(-t^))dtd-d, (97b) 
z J 

S int =/ / ^{^(^^^m^'^OI^-Odtdt'dW, V mt = V + W. (97c) 



where sites are labeled with I and the self-action term V\_(4>i) (|97bj ), that given by Eqs.fl8T|) 
and (91), is separated out. The last term S int describes the two-particle interaction V and 
the effective potential W is caused by averaging over quenched disorder. The potential V 
is assumed to be attractive and takes the standard form ]15| 



Im 

\ E Virnfafa $Mt', #)<i>m(t, #)<t>m(t' , (9* 
Im 

that, in the mean-field approximation, provides the following expression 
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v ~ ~c(t, ^ t, &) x: ^(f, ^)^(*'> - 

~C(t,tf;^#)^(t,tf)^#). (99) 
1 i 

Hereafter v = J2m v im > is the interaction constant, C(t, i?; t', = (<p m (t, i?)</> m (t', $')) 
is the SUSY correlator in the site representation. Averaging over quenched disorder in 
intersite couplings results in the effective attractive interaction || 

W = ~J2 W lrnMt,V)Mt\V')<f>Ut,#)<f>m{t',#'). (100) 
Im 

By analogy with Eq.fl99|) it is supposed that 

in 

w ~ --c(t, t', ■&') x: <f>i(t, ■&'), w = Yl w im > o. (ioi) 

Z I m 

So, the real interaction (|99|) contains both diagonal and non-diagonal in -d and SUSY 
correlators, whereas the quenched disorder averaging results in non-diagonal expression 
( |101|) only. Obviously, within the framework of the replica approach such SUSY structure 
corresponds to the inter-replica overlapping that is responsible for the specific spin-glass 
behaviour [EDI. 



Apart from the above contributions to SUSY action (97) it should be taken into account 



that the quenched disorder in force dispersion results in the additional interaction || 
h 2 

i 

where to is the frequency and the intensity of the quenched disorder 



AS ° = y E / <PU#)Ku)<PU#)foM- (102) 



characterizes the site dispersion of the force fi (overbar denotes the volume average), 
(A<^) 2 = ^5=0 i s mean-squared fluctuation of this force. Then, the mean-field SUSY 
action takes, in the site-frequency representation, the final form 

s = £ / + E / K(#,&)p-dW (104) 

l J In j J III 

with the SUSY Lagrangian 

A(tf) = -</>(•&) {[1 + D(#)] + 27rh 2 5{u) - vS} 4>{$) + (105a) 
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(105b) 



where indexes Z, are suppressed for brevity and the generator D is given by Eq.(|H|). 

In the important case of random heteropolymer the interaction kernels are appeared to be 
of the form fl9~8f ), ( |100| ), but indexes I, m denote wave vectors but site numbers (see 0]). 
So, in this case, the expressions ( P-04j ), (105) can be modified by replacing site indexes by 
wave ones. 



6.2 SUSY Dyson equation 



The Dyson equation for the above SUSY Lagrangian is 

CT 1 =L-£-(v + w)C. (106) 
Here L is defined by Eq.(78) where the first component is 

L = L + vS, L = -(1 + 2nh 2 5(uj)). (107) 



Projecting Eq.( |106| ) along the "axes" (|B"9"1), we come to the key equations written in the 
frequency representation 

G^ 1 + (v + w)G± = L± - E ± (108b) 



where Eq.(]70[) is used. These equations accompanied by Eqs.(96) for the components E, £± 
of the SUSY self-energy function form the closed system of equations for self-consistent 
analysis of non-equilibrium thermodynamic system. 



7 Non— ergodicity and memory effects 



As it is well-known, the memory is characterized by the Edwards-Anderson parameter 
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q = (77(00)7/(0)) (109) 

that being late-time asymptotics of the correlator results in elongation of the structure 
factor: 

S(t)=q + S (t) (110) 
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where the component So(t) — > at t — > oo. By analogy with the elongated structure 
factor ( |110|) , the ergodicity breaking is allowed for by adding the term to the retarded 
Green function 



G?_(w) = A + GLoH. 
The non-ergodicity parameter (irreversible response) in Eq. fllllj ) 

A = Xo ~ X 



'III) 



'1121 



ibo susceptibility Xo = G-(u = 0) and the ther- 
If the latter is defined by the standard formula 
X = 5 (rj) /df ex t with external force f ex t = —f, for determination of the former one has to 
use the correlation techniques discussed in Sect. IV. 



is determined by the adiabatic Cu 
modynamic one x = G-q(lj = 0). 



To do this let us insert the elongated correlators ( |110| ), ( |111| ) into expressions (96). Then 
the renormalized components of the self-energy function take the form 



(113a) 



S±(t) = (^/i 2 + -qj q (A + G ±0 (t)) + E ±0 (t), 
S±o(0 = (fi 2 + X 2 q) S (t)G± (t) + ^S 2 (t)G ±0 (t), 



(113b) 



where S , £±o consist of the terms nonlinear in correlators So, G±o and the terms pro- 
portional to SqA ~ are disregarded. 



In the (^-representation, inserting the Fourier-transform of Eqs. QllO]) , ( |113a|) in the Dyson 
equation ( 108a ), and taking into account Eq. (|107 ) we have 



go 



i 2 1 / 2 I 2 

1 - wxo ~ j [P + Y Qo J qoXo 



u2 2 

h Xo 



So 



;i + s )g + g_ 



1 - [w + (fi 2 + X 2 q/2)q] 



(114) 
(115) 



2 It is convenient to use the unique response function G_ (u) for definition of both susceptibilities 
Xo and \i taking into account that quantities xo = G-(u> = 0), and x = G^{uj — > 0) correspond 



to the equilibrium (macroscopic) and non-equilibrium (microscopic) states. Then Eqs. (p5|) , (|89f) , 
(90), (96), where the correlators should be labeled by index 0, imply the limit u> — > instead of 
the exact equality u = 0. 
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The first of these equations corresponds to 5-terms {uj = 0) that are caused by the memory 
effects, whereas the second one — to non-zero frequencies uj 7^ 0. In the limit uj — > 0, the 
product tends to x 2 , so that the pole of the structure factor ( |115[ ) determines the 

point of ergodicity breaking for the thermodynamic system 



( 2 A 2 \ 



(116) 



Substituting Eq.( |113b] ) into the Dyson equation ( |108b| ) yields the relation for retarded 
Green function 



G-J + 



(v + w)+ (fi 2 + y q j q 



=0 



117) 



where the cu-representation is used. Then, from Eq.( |96cp the equation for the thermody- 
namic susceptibility x = G-(u — > 0) is derived 



1 - x + {v + w)x 2 + —x (x + q) 2 - q 2 



+ 

o 



(x + <?) 3 - <r 



The macroscopic memory parameter go is given by the equation 



118) 



V A 2 \ 2 u2 

Y + % = h, 



119) 



which is obtained from Eqs.(|114l), (|116|) in the limit uj = 0. 



8 Discussion 



Within the framework of the model under consideration the system of Eqs.( |114j ), (|118 



(|116|) , (|119|) , ( |1 1 2| ) provides the complete analytical description of the non-ergodic thermo- 
dynamic system with quenched disorder. Eqs.( |114| ) and ( |118[ ) are similar to the equations 
obtained by Sherrington and Kirkpatrick for determination of isothermal Xo an d thermo- 
dynamic x susceptibilities and corresponding memory parameters go, g M. The equation 
(|116|) defines the point To of ergodicity breaking and Eq.( |112j ) gives non-ergodicity pa- 
rameter A. The above consideration implies that one should distinguish the macroscopic 
quantity g , Xo an d microscopic ones g, x (the former correspond to frequency uj = 0, 
the latter — to the limit uj — > 0). The peculiarity of such a hierarchy is that macroscopic 
values g , Xo depend on the amplitude h of quenched disorder only, whereas microscopic 
ones q, x — on temperature T. Respectively, Eqs. (|114[ ), (|116|) , where the temperature 
T should be taken equal to its value on the ergodicity breaking curve T (h), give the 
macroscopic values go, Xo- What about determination of the microscopic ones g, x-, t ne 
Eq. ( |118j ) must be added by the equation type of Eq. ( |114| ) 



24 



-1 2^/2-^ 2 \2 

1 - w Xo ~ g I ^ + Y q I qXo 



!,2 2 



(120) 



where the memory parameter q is taken as microscopic in character. 



It is well to bear in mind that the field h, anharmonicity parameters A, fi, and interaction 
parameter w, as well as the inverse values of susceptibilities Xo? X have been measured in 
units of temperature T. Further, it is convenient to choose the following measure units: 



Xs 



3\ 3/ V ,3/? 
2J A3' S ~2A 2 '^ 
l ' 2 A 3/i 2 



-1/2 



A. 



2 , q s 



2 A 2 



= u 



121^ 



for quantities T, /i, v, w, x, g respectively. As a result, the key equations take the final 
form: 



1 - uT X ) + (v + w)T X 2 + (X/2T) (T X + qf - q 



+ 



(X/4T) [(T X + g) 3 - g 3 J = 0, 


(122a) 


u;T+(l + g/2)g/2 + /i 2 /g = Xo 2 , 


(122b) 


»T+(1 + g /2) g /2 + /i 2 /g = Xo \ 


(122c) 


[l + (3/4)g ]go + ^T = Xo 2 , 


(122d) 


(l + g )g 2 = 2/i 2 , 


(122e) 


A = uT(xo ~ X)- 


(122f) 



Behaviour of the system is specified by the last parameter u in Eqs.(121), that deter- 
mines the relation between cubic and quartic anharmonicities. In the case u < 1 main 
contribution is due to the quartic term fl81|), whereas at u ^> 1 the cubic anharmonicity 
(|9l| ) dominates. The former limit corresponds to strong quenched disorder h ^> h s , the 
latter — to weak field h <^ih s . 



According to Eq. (|122e|) , the dependence of the macroscopic memory parameter go on 
the quenched disorder amplitude h is governed by the ratio between magnitude h and 
characteristic field h s . The linear dependence 



go = 2 1 / 2 (/i//i), u > 1, 
is realized in the limit h <C h s , whereas the power relation 



(123) 



(124) 
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corresponds to the case of h ^> h s (in Eqs.( |123| ), ( |124|) measured units are restored). The 
dependence qo(h) is depicted in Fig.l. 



For temperatures above the point of ergodicity breaking T the thermodynamic x> Q an d 
adiabatic Xo, Qo values of susceptibilities and memory parameters, as well as Eqs.(122b), 
(122c) coincide, so that Eqs. (122a), (122c) deteremine the dependencies x an d q on tem- 
perature. Accounting Eq.(122d) gives the temperature of the ergodicity breaking T (h) 
as a function of field h (see Fig. 2). The peculiarity of the dependence T (h) is that the 
ergodicity breaking temperature takes non-zeroth value T 00 = Tq(H — 0) at h — 0. Below 
the ergodicity breaking temperature To Eqs. (122a), ( 122b ) give the microscopic values 
of the memory parameter q and the susceptibility x that differs from the macroscopic 
one xo being constant. As a result, we obtaine the typical temperature dependencies of 
susceptibilities x-> Xo as shown in Fig.3. It is seen that, in accordance with Eq.(122a), at 
T < T the thermodynamic susceptibility x 7^ only if temperature is above the freezing 
point Tf which is determined by the condition dx/dT = oo, leading to the equation 



v 



+ w)T f + T fX + q + (3/4) (T fX + qf = 



X 



-2 



(125) 



where Xi 1 are taken at T = Tf. 



The phase diagrams, which depict the ranges of possible thermodynamic states on the 
plane h — T for various values of interaction parameters w and v, are shown in Fig. 2. In 
the limit h = the temperatures of ergodicity breaking and freezing are as follows: 



oo 



v 1 A 2 " 
w < [ H 1 

2w 12 w 2 



1 JL jL^_ 

+ 2w + 12 w 2 



+ 



2w 2 



1/2' 



'126a) 



T f w A{v + w) 1 + 



li 2 + (2/3)A 2 
4(t) + w) 2 



(126b) 



where measured units are restored and the second equality is for /x 2 , A 2 w 2 . So, cubic 
and quartic anharmonicities result in an increase of both ergodicity breaking and freezing 
temperatures. When quenched disorder is large, q , q 2 3> wT , the isothermal susceptibility 
Xo in Eq. (|P22d| ) is small, and Eqs. (122a), (122c), (122d) provide the estimate xo ~ 2/mT . 
Then, for measured quantities one has: 



T w 2 5 / 4 n {h/fi) 1 / 2 , {w/fi) 2 fi < /i< (/i/A) 2 /i; 
T w 2 1 / 2 3 1 / 3 A {h/\) 2 l\ h » WA)V, (w/Xf 2 A 



(127a) 
(127b) 



for jj? ^> A 2 and A 2 ^> /i 2 , respectively. So, the non-ergodicity domain is extended indefi- 
nitely at strong increasing of quenched disorder. As it is shown in Figs.2a,b, and Fig.4a,b, 
the dependencies Tq(/i), Tf{h) is non-monotonous if either w < 0.5 or v > 1. 
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Influence of the interaction parameters w, v and the anharmonicity ratio u on the tem- 
perature dependence of the susceptibility x is illustrated in Fig.5. According to Fig. 5a, 
increasing w causes a decrease of x an d an increase of temperatures T , Tf. The same be- 
haviour is revealed at increasing parameter v (Fig.5b). By contrast, tendency is opposite 
under an increase in u (see Fig. 5c). 

Finally, let us consider behaviour of the non-ergodicity A and memory q parameters that 
are determined by complete system of equations (122). Corresponding dependencies on 
temperature are shown in Fig.6a,b. At the freezing state, where x = 0> the non-ergodicity 
parameter ( |1221j ) linearly depends on temperature because the isothermal susceptibility 
Xo is constant. The appearance of finite value of the thermodynamic susceptibility x above 
the freezing point Tf results in step-like decrease of the value A. With further growth 
of temperature the irreversible response A(T) monotonously decays taking zero value at 
the ergodicity breaking point To (see Fig. 6a). With increasing temperature from to To 
microscopic memory parameter q monotonously decriases, taking minimal value at the 
ergodicity breaking point T . Above this point q{T) increases (see Fig. 6b). It is seen, that 
the quenched disorder encrease extendes the temperature domain of the non-ergodicity 
and causes growth of the memory parameter. In a spirit of generalized picture of phase 
transition it can be attributed to the fact that the microscopic memory parameter q 
above the point T corresponds to a soft mode that transforms to a mode of ergodicity 
restoring below the temperature T . The non-ergodicity parameter A represents the order 
parameter. 

The analytical expressions for dependencies A(T), q(T) can be obtained only near the 
ergodicity breaking curve T (h). For h = and T = T o, from Eq.( [122a| ) assuming that 
< Too — T <ti Too, X ~ Xoo — A/mT 00 , A wXooTw, up to the first order in small 
parameters e = T/T 00 — 1 and A(mxooTdo) -1 we have for measured units 

A = -A e, A,,^{-)' — i— £ <0; (128a) 

» W i + (£ A +5f) (a.) 1/2 

In the case of h 7^ the result for temperature dependence is 

Xo \x^ + >? + qo + 2X0-L0J 

Correspondingly, at the fixed temperature Eq. (|122a| ) gives in the linear approximation 
< g - q < Qo- 

A = B(q-q ), B-^l + ^-t^ + qo + ^XoTo) . (130) 
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A 



For nilpotent representation let us rewrite the Lagrangian ( |10D in the form of Euclidean 
field theory § 

L = K + TT, (A.l) 

where the kinetic k and potential ir energies are 

k = Lpfj — y? 2 / 2 > (A.2) 
7T = (dV/d7])ip. (A.3) 

In order to obtain the nilpotent form ( |13aj ) of the kinetic energy ( |A.2| ), we have to de- 
termine the operator D. The complete form of the dependence of the operator D on the 
nilpotent coordinate $ is presented by the expression 

D = a + b{d/d-&) + cd + d-&{d/d$) , ( A.4) 



where the coefficients a, b, c, d are unknown operators. The substitution of Eqs.(Pl), ([A.4]) 



into Eq. ( |13a| ) and taking into account the properties ( |T2"D leads to the expression (|A.2 ) 



with the following coefficients: 

a = d t , 6=-l, c = 0, d=-2d t , (A.5) 

where c?t = (9/9t is the derivative with respect to time. As a result, the operator (|A.4|) 
takes the form (Q). It has the property 

D 2 = dl (A.6) 



While considering the definitions (0), (|12), (0) it is easy to see that D is a Hermite 
operator. 

Under the infinitesimal transformation 5 = e eD — 1 ~ sD with the parameter e — > 0, the 
values t and i? acquire the additions St = e, S$ = —e which differ in the sign. Considering 
the corresponding field addition \8<f> v \ = elD^Wcp^l, it is convenient to use the matrix form 
for the nilpotent field (|TTD and the operator D: 



v^)' o -a ' dt dt ' 



w = (.:). \ D *\ = I h 9 * = ^- ( A - 7 ) 



According to Eq.( [A.7| ), the change of the order parameter is proportional to a difference 



between the rate of change of the order parameter and the fluctuation amplitude, whereas 
the change of the latter is proportional to its rate with the opposite sign. 
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To prove the equivalence of the term ( |A.3|) in the Lagrangian (|A.1[) and the Grassmann 
potential energy ( |13b|) , let us carry out the formal expansion of the thermodynamic po- 
tential in powers of the component flip of Eq. (|iT|) : 



7T 



5V 

V(v) + T- <P0 
or] 



(A.8) 



Here all the terms of powers higher than 1 are omitted according to the nilpotent condition. 
Using the integration properties (0), we obtain immediately Eq. ( |A.3| ) as it was required. 



In the case of the two-component nilpotent field (p7|), the consideration is fulfilled by anal- 
ogy. For brevity, let us point out the difference between (p7|) and the above-considered 
case (|TT|) only. The corresponding infinitesimal transformation 5 ~ sD, where the gener- 
ator D is given by Eq.(|28]), results in the additions St = 0, 8$ = —e, \6<f>f\ = e\Df\\<pf\ in 
which the matrix form 



v 



-f. 



1 





-d? 



9t =dt 



(A.9) 



is used. The generator Df has the property ( |A.6| ). It is easy seen that the Lagrangians (|10D, 



are invariant under the transformations given by the generators D v , Df, respectively, 
provided that the infinitesimal parameter e is pure imaginary, and the fields rj(r, t), tp(r, t), 
/(r, t) are complex- valued. So, for real fields the two-component representations (|TT|) , (|27|) 
are just convenient approximations. 



The matrices of the transformation between the fields (|A.9 ) and ( |A.7| ) (see Eqs. (|29|) , (^) ) 
take the form 




(A.10) 



Let us consider the four-component SUSY fields (|35|), (|4^). Instead of Eq.( |A.6|) , the 
corresponding couples of operators (57), satisfy the conditions: 



v 2 



0, {V,V} = -2d tl [V,V] 2 = (2d t 



tj > 



(A.H) 



{V v , V f } = {V^, V f } = 0, {Vp, Vj) = -dt, {D f , D v } 



-3d, 



■t, 



where the curly and square brackets denote anticommutator and commutator, respec- 
tively. The generalized anticommutation rules for the operators = X>(±t), = 
D(±t) corresponding to the opposite directions of the time t, read: 



{£>MpW} = {V {±) ,V iT) } = {Vf\vf ] } = 0, 
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{v£\vf} = {vf\vf} 



{V^,vf} = {Vf\vf} = 0, {Vf\vf>} = {Vf\v^} = ±d t . (A.12) 



In Eqs. ([A.ll| ), (|A.12|) the coincident indexes are suppressed. The simplest way to prove 
Eqs. (|A.ll|) , ( |A.12| ) is to introduce the four-rank matrices (see Eqs. (|A.7|) , QA.9Q ): 





f V ^ 


( 


1**1 = 


ill 
















V v J 


V 




( n \ 








, 1^/1 = 








{-fj 


V 



1 0^ 

-2d t 1 



1 0^ 

-$001 



\Vu 



10 
-1 

\0 J 

^001 ^ 
-d t 0-1 

{ 0-dt o J 



(A.13) 



(A. 14) 



The matrices of the transformation between the fields ( A. 14 ) and QA.13 ) take the form 
(cf. Eq-flOq)) 



1 0^ 
10 
10 



y ±dt 1 



(A.15) 



The infinitesimal transformations 5 ~ e~T>, 5 ~ Ds give the following additions: 



^6 = 0, 

5 f e = o, 



-2e6, 



5 V 9 



s v e = o, 



0; 



5 f = E, 6ft = -W, 5 f 9 = E, 8 f e = 0, 6 f t = -6e; 



(A.16) 



\5%\ = el^H^I, \S%\ = \V v \\$ v \e- 
\6Qf\=t\V f \\Q f \, \Uf\ = \Vf\\* f \e. 
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At last, the equation for the SUSY field (j35|) 

/^)^ = ^-?0V> (A.17) 

is obtained by analogy with Eq. (|A.8|) to represent the terms in Eq.(|4]) that contain the 
potential V{i]} in SUSY form. In the case of the field (f44|), the multiplier (p must be 
substituted by — /. 



B 



Following the standard field scheme 0, let us show how the four-component SUSY field 
(pH] ) is split into a couple of chiral two-component Grassmann conjugated fields §±. These 
SUSY fields are obtained from the initial SUSY field under the following transforma- 
tions: 



T± = e ±a , d = eed t , d t = d/dt. 



Accordingly, the generators (^) take the form 
V± = T±V f T T V± = T±V f T T . 



(B.l) 



(B.2) 



Due to the Grassmann nature of the parameter d in the operators T±, it is convenient to 
rewrite (P-2|) in the following form: 



V ± = V f ±[d,Vf], v ± = v f ± 



d,V f 



(B.3) 



where the square brackets denote commutator. In explicit form one has 



V + = d/86 - 26d t , V_ = d/d6; 
V + = d/dO, P_ = 8/ 39 - 26 d t . 



(B.4) 



Apparently, the operators V + , V + coincide with the generators V^, T> v , Eqs. (|37|) . 
According to Eqs. (|4^) , the definitions (B.l|) give 

$± = 77 + 00+^=1=00(17^/), (B.5) 

where the point denotes the derivative with respect to time. The comparison of Eq. ( |B.5| ) 
with the definition (p5|) gives the identity $ + = The action of the operators ( |B.4|) on 
Eq.([B.5|) gets 
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£> ± $ ± =^-0(fi + f) + 299ip, 

-V T <S> T = $ + 9(v- f)-26($, (B.6) 
where the underlined terms concern only the upper indexes of the left-hand parts. 
The chiral SUSY fields are fixed by the gauge conditions || 

= 0, = 0, (B.7) 



which, in accordance with the definitions ( P-4|) signify, that $_ and $ + are independent 
of 9 and 9, respectively On the other hand, taking into account Eqs.( B.6| ), the gauge ( B.7| ) 
results in the equations 

4-Hf + v) = o, 

V> + (t) - f) = (B.8) 

for $_ and $ + , correspondingly Substituting Eqs.( p.8|) into Eq.( p.5| ), the final expressions 
for the chiral SUSY fields are obtained: 

0- = rj + 

<j> + = n + 6ij). (B.9) 



These equations give the non-reducible representations of the SUSY fields (|35|), (|44l) 
under the conditions of the gauge ( |B.7| ). The chiral field <j>+(t) corresponds to the positive 
direction of the time t, whereas 4>-(t) is related to the negative one Q. 



C 

Let us consider the invariance properties of the SUSY action 

S = J[K($(z)) + V($(z)j\dz, K($) = z = {r,t,9,9} (C.l) 

under the Grassmann conjugated transformations 

5$ = ^£ a pWf, 5$ = ^P (a) $e«, (C.2) 

a a 

given by the SUSY generators T>^ , which differ in the time t and Grassmann coor- 
dinates 9, 9. According to Eq. flA.17D , the potential term in Eq. ( p.l| ) is SUSY invariant if 
the kernel V{r]) does not depend on the time t. Up to inessential total time derivatives, 
the Grassmann conjugated variations of the remaining kinetic term 
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(C.3) 



can be rewritten in the form 



SK = -J2V (a) [(£>$)(£>$) 

2 -T, 



(C.4) 



provided that the anticommutators {£>,£>(")}, {£>,£> {q) }, {£>,£> (q) } are either 

equal to zero or proportional to the derivative with respect to the time d t . According to 
Eqs. flATTD , (|A~T^) such condit ions are fulfilled if only the SUSY generators D", T> 
either coincide with the initial operators V, V, or are reduced to the transformed op- 



erators P±, V± determined by equations of (|B.2|) type, or correspond to the opposite 

time directions T>^\ T>^±\ Being reduced to the derivatives with respect to the time t 
and Grassmann coordinate 9, 8, these operators inserted into Eqs.( p.4] ) give, as it was 
required, zero for the variations of the corresponding action ( |C.1| ). 



Among the above-mentioned generators, the following ones 



V 



(-) 



d_ 



d 

de + 29 dl 



(C.5) 



and their anticommutator {v[_ \ v[_ = 2d t (see Eqs.( A.12 )) are of a special interest 
for us. The operators ( |C.5| ) are a result of the double action of the transformation T_ 
on the initial generators T)^, Eqs. (|3~7D , that gives the generators D_, £>_, Eqs.( |B.3| ), 
corresponding to the opposite time directions. Therefore, the generators P_ ■* = X>_(— £), 

given by Eqs.( |C.5| ) are related to the initial ones T> v {t), D^it) and play a 
significant role hereinafter. 

Due to the standard manner, it is easy to show that the above conditions 5S = 0, 5S = 
give rise to the Ward identities M 



$>J a) r(«)({^}) = o, y^rM^}) = o 



(C.6) 



i=l 



i=l 



for SUSY n-point r^ n ^ proper vertices type of the 2-point supercorrelator C(z2, Zi), 
Eq.([53]) and the self-energy superfunction T,(z 2 , Z\). Obviously, under T>^ = d t , T>^ = 
^ conditions ( |C.6| ) mean that above-mentioned supercorrelator depends on t 2 — ti and 
62 — 0i differences only: 



C v (z2,zi) = Sfa-h) + 



G-(t 2 -ti)6i 



(C.7) 



33 



where the space dependence is suppressed, for brevity. In accordance with the SUSY field 
definition (^), one has: 

s{t 2 - h) = {vihHti)) , 

G+(t 2 - h) = (ipfoWh)) - t 2 ) = (^(t 2 )V(*i)) fifa - h), (C.8) 

G-(t 2 - h) = (vihMh)) #(t 2 - h) = (^(t^ih)) #(t 2 - h) 

where the step function = 1 for t > and = for t < 0. By virtue of the causality 
principle, the advanced Green function G + (t 2 — 1\) which is the factor before 8i6 2 , vanishes 
at ti < t 2 , as required. On the other hand, the retarded Green function G_(t 2 — ti) = 
at t\ > t 2 , to be the coefficient of 6 2 6i- Moreover, the symmetry condition C(z 2 ,Zi) = 
C(zi,z 2 ) gives rise to the equations S(t 2 — ti) = S(ti — t 2 ), G-{t 2 — t\) = G + (t\ — t 2 ). 
Inserting the operator T)[_ ^ into the Ward identity (|C.6|) results in the equation 



2S(t) = G + (t) - G_(t) (C.9) 

that is the fluctuation-dissipation relation. All the above statements hold for the self- 
energy function £(z 2 , Z\) with the components T>(t 2 —ti), £±(t 2 — ^i) which replace S(t 2 — 
ti), G±(t 2 — ti), correspondingly. 



It is worthwhile to point out specially the relations in Eqs. (|G.8| ) which connect the Bose 



correlators for the components rj, (p, and the Fermi ones for components ip, ip. The above- 
used Ward identities ( |G.6| ) allow to obtain these relations as a trivial consequence of 
the SUSY field definition ([351). But such equations can be obtained also more simply. 
Indeed, the Fermi correlator (ipip^ is equal to {{5 2 V / 5r] 2 )~ l ) in accordance with Eqs.(f21~D, 
(|34]). On the other hand, using the susceptibility definition and Eqs. (|22|) , fl26|) we have 
{r)(p) = {Sr)/S(p) = {(S^/St])^ 1 ) = {(S 2 V / 'St] 2 )" 1 ) for the Bose correlator Q.E.D. 

Obviously, to pass to the correlator of the two-component field (pTTj) , it is necessary to 



replace the factors (6 2 — 6i)6 2 , (6 2 — 6i)6i in Eq. (|C.7|) by the nilpotent coordinates $ 2 , 
— respectively, and to omit the term with i? 2 #i- The SUSY correlator Cf (z 2 ,z{) = 
T^(z 2 )T^(zi)C ip (z 2 , Z\) corresponding to the superfields (|27|) and (£0]), takes the form 



C f (z 2 , zi) = S + 9 2 9 2 m + + e x 6 x m- - 6 2 e 1 G- - 9 1 6 2 G, 



(CIO) 



where the arguments t 2 — ti are suppressed in the factors S, m±, G±. Moreover, in view 
of Eq. fl22|) , new functions are introduced (cf. Eqs. (|G.8| )) 



m+(t) = ( V (-t))^(-t)U, m_(t) = {r } (t))0(t)f eKt , f ext = -f (C.ll) 
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to represent the connection between the averaged order parameter (i](t)) and the external 
force / ex t = — / (note that the latter is switched at time t = and remains constant). 
The correlators ( |C.11| ) are related to the Green functions as follows: 

G±{t) =m±{t)±S{t) (C.12) 
and possess the symmetry condition m+(—t) = m_(t). 
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Captures 



Fig.l Dependence of the macroscopic memory parameter q on the quenched disorder 
intensity h. 

Fig. 2 Dependencies of the ergodicity breaking temperature T (solid line), and the freezing 
temperature 7/ (thin line) on the quenched disorder intensity h (u — 0.5, v — 0) for 
different values of the effective interaction parameter: a) w=0.5; b) w—0.2. 

Fig. 3 Temperature dependencies of the thermodynamic and adiabatic susceptibilities x 
and Xo f° r: a ) different values of the quenched disorder intensity h (curves 1, 2 correspond 
to h — 0,4) at w — 0.5, u = 0.5, v — 0; b) different values of the effective interaction 
parameter (curves 1, 2 correspond to w — 0.5, 0.2) at h — 4, u — 0.5, v — 0. 

Fig. 4 Dependencies of the ergodicity breaking temperature T on the quenched disor- 
der intensity h and: a) effective interaction parameter w at u — 0.5, i> = 0; b) proper 
interaction parameter v at u — 1, w — 0.5. 

Fig. 5 The shift of the temperature dependencies of the thermodynamic susceptibility xCO 
caused by variation of: a) effective interaction parameter w at u = 0.5, v — (curves 1, 
2, 3 correspond to u> =0.5, 1, 1.5); b) proper interaction parameter v at u = 0.5, w = 0.5 
(curves 1, 2, 3 correspond to t>=0, 1, 1.5); c) anharmonicity ratio u at w — 0.5, v = 
(curves 1, 2, 3 correspond to w =0.5, 1, 1.5). 

Fig. 6 Temperature dependencies (u — w — 0.5, v — 0) of: a) non-ergodicity parameter 
A; b) microscopic memory parameter q. (curves 1, 2 correspond to h =0, 4). 
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Fig.l 




Fig.2 




Fig. 6 



